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Abstract 



I identify a source of ^-qcd/^q corrections to the assumption of quark- 
hadron duality in the application of heavy quark methods to inclusive heavy 
quark semileptonic decays Q — * qivi- These corrections could substantially 
affect the accuracy of such methods in practical applications and in partic- 
ular compromise their utility for the extraction of the Cabibbo-Kobayashi- 
Maskawa matrix element V c b- 
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Although the classic application of heavy quark symmetry is in the exclusive semileptonic 
decays of heavy quarks [[[]] , there has also been substantial work on using heavy quark effec- 
tive theory (HQET) [[J to systematically improve decay predictions for inclusive semileptonic 
decay rates of heavy hadrons induced by an underlying Q — > qti>i decay [0-0] • In these in- 
clusive applications, decays are treated in an operator product expansion (OPE) which leads 
via HQET to a I/uiq expansion in which the leading term is free quark decay, 1/mq terms 
appear to be absent, and terms of order l/rrin take into account such effects as the difference 
A = tub — mi, between m# (the mass of a B meson) and my, (the mass of a b quark), the 
nonzero kinetic energy (— Ai/2ra&) of the heavy quark, small residual spin-dependent inter- 
actions (2X 2 / rrif) = rriB* — mg) of the heavy quark, etcetera. Although these calculations 
have become very sophisticated HJ5|], it is widely appreciated ||[7j that there remains a basic 
unproved hypothesis in their derivation: the assumption of quark-hadron duality. It is the 
accuracy of this assumption that I want to call into question here. 

While supposedly valid for any semileptonic decay Q — > qtvg of a heavy quark Q, re- 
cent applications have centered around the hope that this approach offers an alternative to 
the classic exclusive methods for determining V c b, and I will accordingly focus most of my 
remarks on the case b — > dv^ where both quarks are heavy. In inclusive b — > cli>i decays, 
which materialize as B — > X c £u£, about 65% of the X c spectrum is known to be due to the 
very narrow ground states D and D* . The relatively narrow sj e = | + states || 1^2(2460) 
and _Di(2420) account for perhaps another 5% of the rate, and it may be assumed that 
the remaining rate involves decays to higher mass resonances (quarkonia and hydrids) and 
continua ||. The HQET-based inclusive calculations predict continuous X c spectra which 
are assumed to be dual to the true hadronic spectrum (see Fig. 1). 

A picture like Fig. 1 might lead one to dismiss the duality approximation since the inclu- 
sive spectrum clearly does not meet the usual requirement that it be far above the resonance 
region [|K]]. I.e., normally the accuracy of quark-hadron duality would be determined by a 
parameter Aqcd/E where the relevent energy scale E is the mean hadronic excitation en- 



ergy Amx c = m-Xc ~ m D- However, as first explained by Shifman and Voloshin JTTJJT^] , this 



is not the expansion variable in this case: duality for heavy-to-heavy semileptonic decays 
sets in at threshold since even as 5m = wi;, — m c (and therefore AmjJ approaches zero, as 
rrib — ► oo the heavy recoiling c quark has an energy much greater than Aqcd so that its 
final state interactions can be ignored. In the small velocity (SV) limit, it must therefore 
hadronize with unit probability (up to potential h.QCD/ m Q corrections) as D and D*. This 
"cannonball" approximation is in fact an essential part of the physical basis of the HQET 
expansion in I/uiq. Thus the issue is not whether duality holds in semileptonic heavy quark 
decays, but rather how accurately it holds. 
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Fig. 1: A sketch for b —> c semileptonic decay of the continuous inclusive recoil spectrum 
of the OPE calculations (smooth curve) compared to the known hadronic spectrum (shown 
as individual resonance lines). 

Let me begin my discussion of the accuracy of inclusive methods with an overview. Up to 
caveats regarding the unknown accuracy of the assumption of duality, the OPE indicates that 
inclusive calculations should be accurate up to corrections of order A.q CD / toq. Here I will 
identify a source of duality- violation which leads to Aqcd/itt-q corrections. It is revealed by 
considerations of a Bjorken sum rule which may be viewed as an extension of Shifman- 
Voloshin duality to arbitrary recoils. Bjorken's sum rule guarantees that, as m& — > oo, 
duality will be enforced locally in the semileptonic decay Dalitz plot of rate versus w — 1 and 



Ei (where w = v ■ v' is the usual heavy quark double- velocity variable and Ee is the lepton 
energy). For regions of the Dalitz plot for which w — 1 is not large (and in b — > c decay 
nearly the whole Dalitz plot satisfies this condition), the Bjorken sum rule explicitly relates 
the loss of total rate from the "elastic" s# e 
to the turn-on of the production of sj e = \ + and | + states [H 
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Fig. 2: The exact compensation by inelastic channels of the fall of the elastic rate in the 
linear region as m;, — > oo. 



In particular, in this region the Isgur-Wise function may be taken to be linear: 
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In that case, for fixed r = m c /mfe, as — > oo inelastic sj e = | + and | + channels open up 
to give a semileptonic rate that would exactly and locally compensate in the Dalitz plot the 
loss of rate from the elastic channels due to A, m . I.e., if 



then pip 



iyn' 



J2~pinclusive 
quark 

dwdEe 



K(w, Ei) 



(2) 



M'pinclusive 
u hadron 



K{w, E £ 



w + 1 



dwdEe 

as mb oo and for p 2 (w — 1 



Eki m) (l)| 2 + 2£|rf(l 

2 



2 

max = p 2 - ~2 r ' << 1- Since according to (1) 



Cp)z 

2 



(3) 



(^)lCH| 2 ^l-2pU^- 1 ) > ( 4 ) 

we see explicitly the compensation of the dynamical part of the fall of the elastic channels 
by the onset of inelastic states (of all sorts) with sj e = | + and | + . This situation is sketched 
in Figure 2; if it were applicable to b — > c decays, then quark-hadron duality would be exact. 

Having established conditions for its validity as m& — > oo, it is easy to see why one should 
be concerned about quark-hadron duality for b — > c decays. For fixed r, w — 1 lies in the 
fixed range from to (1 — r) 2 /2r, and as m& — > oo any given hadronic threshold collapses 
to the point w — 1. However, for finite m& there is a gap in w — 1 in which the rate to the 
elastic \ channels falls by Kqcd / m Q terms but the potentially compensating excited state 
channels | + and | + are not yet kinematically allowed. More precisely, if mo** is the mass 
of a generic charmed inelastic state, then t** = (m B — m D **) 2 would be the threshold in t 
for this state, corresponding to a value of w — 1 in the quark-decay Dalitz plot of 

t — t** A 

m tm ~(l- r )— (5) 
2m\ ) m c m c 

where t m = (m B — m D ) 2 ~ (m b — m c ) 2 and A = m D ** — m D . Since A ~ 500 MeV 
and (w — I) max — 0.6, this region covers more than one third of the Dalitz plot and the 
compensation is very substantially delayed: see Figure 3. Eqs. (5) and (1) show that this 
effect is of order Aqcd/it^q, seemingly at odds with the OPE result. 

Despite this apparent contradiction, there is actually no inconsistency: the OPE result 
that the leading corrections to the inclusive rate are of order k 2 Q CD /m 2 Q can still be valid as 
derived in the limit of large energy release in the b — > c transition, while Aqcd/^q effects 
can arise for energy releases of the order of Aqcd due to a finite radius of convergence of 
the OPE. The main purpose of this paper is indeed to call attention to this effect, and to 
use the quark model to estimate its importance. 
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Fig. 3: An overlay of the Dalitz plots for B -> £> ( *W e , B -> L>**(2500)ez/ e , and 5 -> 
.D**(3000)ez/ e . The -D**- 1 mass is taken as the hyperfine average of the D and D* masses; the 
two D** masses are chosen for illustrative purposes. 

The basic issues can be most easily exposed by considering [Q spinless quarks coupled to 
a scalar field <fi of mass fi, and by studying the decay b — > c<p with weak coupling constant 
g. Differential semileptonic decay rates have a more complex spin structure, but otherwise 
correspond to the case ll = \ft; total semileptonic rates correspond to a weighted average 
over kinematically allowed li but, as we shall see below, this averaging does not change the 
essentials of the problem. In our simplified case 

r(6 - o0) = p% (6) 

where pji = [(m, — raf) 2 — /i 2 ] 1//2 [(mj + mj) 2 — /i 2 ] 1 / 2 /2m i is the momentum of from the 
two-body decay of mass rrii into masses rrif and \i. 

To compare Eq. (6) with a hadronic world, we can use the quark model framework of 
Refs. irSIJTBIl where 



,2„ 



T(B -> = ^^(^^)|^(^) B )| 2 (7) 
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with 



£ (n) (v D( n) B ) = J d 3 p(p* D(n) (p)(f) B (p-m d v D(n)B ) 



and where 



(8) 



V D (n) B = P D (n) B /m D (n) 

is the recoil velocity of D^ n \ These results may be used to estimate 

r(6 -> C 0) 

Using simple harmonic oscillator wavefunctions in (8) gives 
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where p B . B = + &d)-> £db is the B ^ D form factor, and is the form factor for 

transitions into the lowest £ = 1 excited state with me = i. If we define an order Aqcd/ttiq 
expansion parameter 



_ md _ m d (m b - m c ) 
fi_ m h m c 
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and the "scaled energy release" 



nib — m c — fi 



(14) 



where A = mo** — mo, then for small T* (i.e., small velocities) 



\i DB \ 2 = l-eT* + 0{e 2 ) 



(15) 



and, defining \£ D ** B \ 2 = £i |^« B | 2 > 



e(T* - 1) + 0(e 2 



(16) 
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Since in this limit 
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we can obtain a model H for R by truncating the sum over n after the first D**\ 
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,D** — 
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(T* - l) 3 / 2 



0(T* - 1) 
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wherein we have made the weak binding approximation (3 « rrid and shown explicitly the 
two thresholds at T* = and T* = 1. As T* — > oo, the complete expression for i? will 
display a tower of states (each with their appropriate threshold factors) being produced 
with a strength given by a power series in T* with coefficients which are power series in 
generalized e-coefficients of order KqcdI^q- For duality to be valid, each term T* n for 
n > must have zero coefficient. The OPE further states that the constant term T*° is of 
the form 1 + 0(A 2 QCD /m 2 Q ). 

While extreme, this truncation of the complete expression for R has the properties that: 

1. At T* -> oo, it is of the form 1 + 0(e 2 ) + 0(e/T*) as required by the OPE. 

2. There are no other terms of order 1, e, or eT* possible beyond those shown: a more 
accurate treatment of T(B — ► D<p) could only generate e 2 , e 2 T*, e 2 T* 2 , ... terms; a more 
accurate treatment of T(B — > D**<p) could only generate e 2 T*, e 2 T* 2 , ... terms; and all 
higher states first make a contribution at order e 2 T* 2 or higher. Conversely, we note that if, 
for example, e 2 T* 2 terms are retained, they must all cancel exactly or the requirements of 
the OPE would be violated as T* — > oo. 

3. As Am = m b -m c ^ 0, R? ** -> 1 + 0( AQCD . Am ) as required pf . 

m i, — 

4. Near T^ aa , = m& — m c , eT^ aa , is in general large. This observation corresponds in the 
usual language of heavy quark symmetry to the statement that the natural scale of the slope 
p 2 of the Isgur-Wise function is of order unity. It is also consistent with the experimental 
observation that |£d_b| 2 has dropped to less than half its value between w — 1 = and its 
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maximum value. Given this, the extension of Eq. (19) to higher orders in T* will require a 
"conspiracy" of the entire spectrum of possible hadronic final states. Nevertheless, we can 
use Eq. (19) across the full range of T* as an indicator of the Aqcd/ttiq effects arising from 
the order 1 and order T* terms in the expansion of R; this corresponds to a "best case" 
assumption that duality is locally perfect for the terms T* n with n > 1. 

While this simple example clearly demonstrates the existence of the claimed duality- 
violating Kqcd/^q effects for finite T*, it remains to discuss their quantitative importance. 
Before doing so, I will introduce a number of simple variants of this prototypical model. 
The first corrects an idiosyncracy of the simple harmonic oscillator model: in it the first D** 
state saturates the Bjorken sum rule, i.e., completely compensates the — eT* elastic term. 
In contrast, in the ISGW model where harmonic oscillator wavefunctions are simply used 
as a variational basis, {3d** ^ Pd and this saturation does not occur. This is typical of the 
general case where (in the narrow resonance approximation) the total resonant P-wave term 
is of the form 

T^Y^UT* -tlf' 2 9{T* -e n ) (20) 

J- n 

with J2n fn — 1 and t* being the threshold for channel n. As T* — > oo, these contributions 
automatically cancel the —eT* term from the elastic form factor, and constrain the 0(e) 
correction: 

RSI... = [1 + let* - eT*]9(T*) + ' £ f n (T* - t* n f 2 9(T* - t* n ) , (21) 

z J- n 

where 

i* = E fnt* (22) 

n 

is the weighted average threshold position. Note that since some T* exceed T^ ax , R±+ 2 +... 
cannot heal to unity in the physical decay region. 

As described above, both i?f* and Ri" 2+ are "best case" truncations which assume 
exact cancellations of e 2 T*, e 2 T* 2 , ... terms. While sufficient for the purposes of this study, 
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I note that it is straightforward to recursively "construct duality" to the required order in 
e to any finite order in T* . For example, for a simple harmonic oscillator spectrum one can 
easily construct 

exp{-eT*) / 3 



R ho = % /2 ] ([i + Ut*^ 9 (t*) 



1 2 16 32 1024 JV ; v ; 

+le 2 [l + 7 -e + ^e 2 + ^e 3 + ...](T* - 2fl 2 6{T* - 2) 
21 2 A lb 

+^e 3 [l + 9 -e + ^e 2 + ...](T* - 3)^9(T* - 3) 

+ ^e 4 [l + ye + ...](T* - 4) 9 / 2 #(T* - 4) 

+ ie 5 [l + ...](T* - 5) n / 2 ^(T* - 5) + ...) , 
5! ' 7 

where the ellipses denote terms of order e 6 T* n with 1 < n < 5 and all terms of order 
e m T* m and higher with m > 5. This truncated expansion is accurate even at T* = 5 up to 
corrections of order Aq CO /jtiq; as we will see below, higher values of T* are probably not 
physically relevant. 

The models just introduced are all based on the duality of b — > c<p to a tower of cd res- 
onances. While the thresholds associated with such towers are a source of duality-violating 
^■qcd/itT'Q corrections which must be a cause for concern in comparing inclusive calculations 
with experiment, I am even more concerned about processes which could give a nonpertur- 
bative high-mass tail to the recoil mass distribution. The convergence to unity of R in the 
former case would be controlled by an expansion in l/T* ~ A/(m& — m c — fi), i.e., T must 
be large compared to the single resonance scale A. Since A/T max ~ 1/6, there are reasons 
to be cautious, but perhaps not alarmed. However, I believe that there is another "harder" 
effective scale in low energy hadron structure which I identify with the constituent quark 
size r q . While, like A, l/r q must be "of order Aqcd \ empirically it is much larger: the 
constituent quark model makes sense only if the quarks are small relative to hadronic radii. 



Quantitative estimates from spectroscopy [18| indicate that this scale is indeed l/r g ~ 2 



GeV. This is in turn potentially very dangerous for duality when T max is only of order 3 
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GeV as in b — > c decays. 

Let me mention two concrete examples of how this scale could lead to a high-mass tail 
to the recoil mass distribution which would delay the compensation required for duality. 
If the glue around a constituent quark is indeed very compact, then high recoil momenta 
of the c quark are required if it is to be left behind, i.e., before the current c quark can 
carry away all of the energy of the underlying b — > c0 transition. This effect would lead to 
a convergence problem for the resonance contributions to R associated with the excitation 
of hybrid mesons. The second example is probably more dangerous: such a scale could 



lead to the nonperturbative production of high-mass nonresonant states [19[], e.g., an X C Y 
mass spectrum for B — > X c Y(j) extending from X C Y threshold up to masses of order nix c + 
my + l/r q ~ 5 GeV ~mg. We should therefore be concerned that a substantial fraction 
of nonresonant production is unavailable to participate in the compensation required for 
duality to be realized in b — > c decays. 

For a crude estimate of the effects of a high-mass hybrid contribution, I take a simple two- 
component resonance model consisting of "normal" cd resonances with t* d and cd hybrids 
with i* hybrid substantially larger. If we assume that the latter are responsible for a fraction k 
of p 2 dyn , then we would have 

R hybrid = ^ + _ eT *]0( T *) 



W Jtf e(T* - 1» ) , (23) 



with t* — (1 — K)i* cd -+ ^hybrid- Since both experimentally and theoretically the hybrid mass 
spectrum begins about 1.5 GeV above the ground state, this provides a minimum value 
for i* hybrid . Since their "hard" production mechanism will raise their mean T* above this 
minimum, the effects of their postponed onset could be serious. 

As already implied, I believe that the postponed nonresonant contributions are an even 
more serious cause for concern. While the model (23) for the hybrids could also be used as 
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a template for a crude model for nonresonant states (their contribution will in first order be 
controlled by a parameter t* nr which would replace i* hybrid in Eq. (23)), a more appropriate 
model would be 



R nr = [1 + -et* - eT*]6(T*) 

(T* — f*_") 3 / 2 



rptfl/2 



+Ae / dt*p(t* 



/rp* _ £*N3/2 
ji*l/2 



(24) 



where A is the fraction of p 2 dyn due to nonresonant states and p(t*) is the appropriate nor- 
malized spectral function (J*™ dt*p(t*) = 1) which begins at T^ in but drops off very slowly 

min 

with a scale determined by l/r„. In this situation, t* — (1 — A)t*7 + A J™ dt*p(t*)t*. 
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Fig. 4: Four resonance models of the approach to duality: (a) * (with the baseline value 
e = 1/6), (b) i?f (with e = 1/3 corresponding to T^ a2 , = 3), (c) Ri" 2+ (with t* n = n and 
/ n = (|) n so that t* = 2), and (d) 

I now turn to rough quantitative estimates of duality- violating effects using the models I 
have presented. Since experimentally p 2 ~ 1, (w — l) max ~ 1/2 and — m c )/A ~ 6, I will 
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take e = 1/6 and T^ ax = 6 as a realistic baseline for all of the following examples. However, 
in keeping with the main message of this paper that inclusive results must be interpreted 
cautiously until duality-violating effects are better understood, I will choose pessimistic 
values for other parameters of the models and for variations around these baseline values. 

Figure 4 shows four different examples of resonance compensation: R® (with both the 
baseline value e = 1/6 and for e = 1/3 corresponding to T^ ax = 3, i.e., to using A e ff = 2 A 
for the mean location of the P-wave strength), R^l + (with t* n = n and f n = (|) n so that 
t* — 2), and R . These examples show that duality- violation of the order of 10% could 
easily arise in b — > c decays from the delayed onset of resonances with a scale A ~ 500 MeV. 




Fig. 5: Four examples of the effects of a nonperturbative high mass tail on the recoil mass 
spectrum: (a) R h v bnd (with k = 1/10, t*^ = 2, and i* hybrid = 4), (b) As in (a), but with 
^hybrid = 6, (c) R nr (with A = 1/5, t* d - = 2, T^ in = 2, and s = 4), and (d) As in (c) but with 
s = 6. 
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Figure 5 shows the potentially even more dangerous effects of there being a second scale 
in this problem much larger than A. Here I once again give four examples: R h v brtd (with 
k = 1/10, i* c3 = 2, and t* hyhrid = 4 and 6) and R nr (with A = 1/5, i* cg = 2, and 

p{e) = lexp{ T ^~ t * )9{t* - T* mm ) (25) 
s s 

with T^ in = 2 and s = 4 and 6). Note that the choices made for t* hybrid and s are based 
on the hypothetical scale l/r q ~ 2 GeV. The choices k = 1/10 and A = 1/5 are based on 
estimates §| of the strengths of these contributions, but values of this order are certainly 
reasonable (e.g., A is a 1/N C effect). These examples show that it is difficult to avoid the 
conclusion that nonperturbative high mass contributions to p 2 could lead to substantial 
duality- violating Aqcd/i^q effects in b — > c decays. 

Although our main focus has been on heavy-to-heavy transitions, the physics issues 
raised here (if not their explicit forms) are also relevant for Q — > qtvi transitions. Before 
concluding, let me therefore point out a simple application of the OPE to inclusive heavy-to- 
light transitions where it seems certain to me that they will fail: Cabibbo-forbidden charm 
decays. (Even though such decays might be an unimportant application of the inclusive 
calculations in practice, they provide a valid theoretical testing ground for their accuracy.) 
In particular, consider the c — > diu£ decays of the D° and D + . They will be dominated by 
the channels D° — > -n~ivi and p~ivt and by D + — > tt ^, rjlvi, r(tvi, p lvz, and ulvt,. Since 
the OPE corrections in the D° and D + are identical, their Cabibbo-forbidden semileptonic 
partial widths and spectral distributions are predicted to be identical. However, simple 
isospin symmetry implies that T(D + — > n^Ivg) = ^T(D° — > ir'ivi), so the inclusive Cabibbo- 
forbidden rates can only be equal if T(D + — > rjivi) + T(D + — > rfive) = T(D + — > 7r°£z/f). 
In many models this latter relation would be true if m v = = m n , since it is rather 
natural for the squares of matrix elements to satisfy its analogue. However, with real phase 
space factors, this relation is typically badly broken. Since Cabibbo-forbidden decays, like 
their Cabibbo-allowed counterparts, will receive little P-wave compensation, I expect this 
prediction to fail. 
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Finally, I note that the duality- violating effects I have highlighted here will have an effect 
on the long-standing B semileptonic branching ratio puzzle [ptfl . Since the hadronic mass 
distribution in b — > cud is weighted toward higher masses than the leptonic mass distribution 
in b — > c£p£, the ratio of these two rates will be changed. 

In summary, I have shown here that hadronic thresholds lead to AgcD/ m Q violations of 
duality in b — > c decays which do not explicitly appear in the operator product expansion. 
Since such violations cannot appear as the b — > c energy release T — > oo, there are "con- 
spiracies" (i.e., sum rules) which relate hadronic thresholds and transition form factors. As 
emphasized by Bigi, Uraltsev, Shifman, Vainshtein, and others fU||7]] , these relations tend to 
compensate the otherwise extremely large Aqcd/ m Q effects even at small T. In this paper 
I have displayed several models of such hadronic compensation mechanisms which indicate 
that these duality- violating h.QCD/ m Q effects could nevertheless be very substantial. While 
the examples I have selected are perhaps pessimistic, they indicate that these effects must be 
better understood before inclusive methods can be applied with confidence to heavy quark 
semileptonic decays. 
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